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Abstract: Martin has proposed a scenario dubbed "compressed supersymmetry" (SUSY) 
where the MSSM is the effective field theory between energy scales M„eak and Mqut, but 
with the GUT scale 5J7(3) gaugino mass M3 <S Mi or M2. As a result, squark and 
gluino masses are suppressed relative to slepton, chargino and neutralino masses, leading 
to a compressed sparticle mass spectrum, and where the dark matter relic density in the 
early universe may be dominantly governed by neutralino annihilation into ti pairs via 
exchange of a light top squark. We explore the dark matter and collider signals expected 
from compressed SUSY for two distinct model lines with differing assumptions about GUT 
scale gaugino mass parameters. For dark matter signals, the compressed squark spectrum 
leads to an enhancement in direct detection rates compared to models with unified gaugino 
masses. Meanwhile, neutralino halo annihilation rates to gamma rays and anti-matter are 
also enhanced relative to related scenarios with unified gaugino masses but, depending 
on the halo dark matter distribution, may yet be below the sensitivity of indirect searches 
underway. In the case of collider signals, we compare the rates for the potentially dominant 
decay modes of the ti which may be expected to be produced in cascade decay chains at 
the LHG: ti — > cZi and ti — > bWZi. We examine the extent to which multilepton signal 
rates are reduced when the two-body decay mode dominates. For the model lines that 
we examine here, the multi-lepton signals, though reduced, still remain observable at the 
LHC. 
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1. Introduction 



Models of particle physics with weak scale softly broken supersymmetry are well-motivated 
by both theory and experiment. On the theory side, they stabilize the scalar sector that 
plays an essential role in the spontaneous breaking of electroweak symmetry, allowing a 
sensible extrapolation of particle interactions over many orders of magnitude in energy. 
On the experiment side, supersymmetric models naturally accommodate, i) gauge cou- 
pling unification, ii) a mechanism for electroweak symmetry breaking due to a large top 
quark mass, in) a light Higgs scalar and decoupled superpartners in accord with preci- 
sion electroweak measurements and iv) a neutral weakly interacting particle that 
thermal Big Bang relic, account for the observed cold dark matter (CDM) in the Universe. 

In spite of the accolades, supersymmetric theories suffer from new problems not present 
in the Standard Model (SM). There are the big issues such as the flavor and the CP 
problems, as well as the fact that baryon and lepton numbers can potentially be violated 
at large rates. We have nothing new to say about these, and will evade these in the usual 
ways. A much less serious objection is the "supersymmetric little hierarchy problem" which 
simply states that the value of the parameter —m'j^^ (renormalized at the TeV scale) can 
be ~ only if there are cancellations at the percent level, once experimental constraints 
on sparticle and MSSM Higgs scalar masses are incorporated. Another potential problem 
is that in many supersymmetric models, the lightest SUSY particle, usually the lightest 
neutralino, is bino-like, with a typical thermal relic density considerably larger than the 
measured CDM density ^cDMh'^ ~ 0.1 for sparticle masses larger than ~ 100 GeV. 

Recently, Martin has observed that the latter two issues are ameliorated in a sce- 
nario that he calls "compressed supersymmetry". Within this framework, it is as- 
sumed that the MSSM is the effective field theory between M^eak and Mgut- As in 
the mSUGRA model, universal scalar mass parameters are adopted at Q = Mgut but 
non-universal gaugino mass parameters are allowed. Specifically, Martin notes that if 
3M3(GUT) ~ M2(GUT) ~ Mi (GUT), the fine-tuning required to obtain small values of 
|m^^ | is considerably reduced. The low value of M3 results in a SUSY spectrum where 
physical squark and gluino masses are closer in mass to uncolored sparticles than in models 
such as mSUGRA with unified gaugino masses, where one expects rriq ~ nig » ''^Wi' '^^^s 
the SUSY spectrum is "compressed" relative to models with gaugino mass unification. 

Of particular interest to us are solutions with a compressed spectrum where the top 
squark ti is particularly light. In this case, if the neutralino annihilation channel ZiZi tt 
is kinematically accessible in the early Universe, its reaction rate suffers no propagator 
suppression because of the light t— and u— channel stop exchange, and can lead to a 
neutralino relic abundance in accord with WMAP, even though the neutralino remains 
largely bino-like. In addition, as noted above, the low third generation squark masses feed 
into the evolution of the soft SUSY breaking Higgs mass rn'jj^ , causing it to evolve to much 
smaller (in magnitude) negative values than in the case of unified gaugino masses. Since 
— m^^(weak) ~ /i^ the little hierarchy problem is less severe than in models with unified 
gaugino masses. 



- 1 - 



Martin has shown that the compressed SUSY scenario is vahd provided that 



+ 25 GeV ~ ^ r 



mt + 100 GeV, 



(1.1) 



mj +100 GeV 

^1 



(1.2) 



where the lower hmits above are imposed so that annihilation of neutralinos into top pairs 
is allowed at rest, and to reduce the impact of ti-Zi co-annihilation, while the upper limits 
should be viewed as soft. He displays an explicit case where the GUT scale gaugino masses 
are related according to 



which can occur in models where the SUSY breaking F-term that seeds SUSY breaking 
gaugino masses transforms as a linear combination of a singlet and an adjoint field of the 
unifying SU{5) group. The trilinear soft SUSY breaking term Aq is set either to —Mi or 
— 0.75Mi. Since the ii — Zi mass gap is small in compressed SUSY, Martin recognized 
that two cases emerge which are relevant to LHC searches: one is characterized by when 
ti — > cZi is the dominant top squark decay channel, while the other has a large enough 
mass gap that ti bWZi can compete, and perhaps dominate, the two-body decay. 

In fact, this whole scenario appears closely related to scenarios first pointed out by 
Belanger et al.'^ and independently by Mambrini and NezriQ and subsequently exam- 
ined in detail in Ref. where a reduced GUT scale gaugino mass M3 leads to a small 
H parameter, and ultimately to a mixed higgsino-bino Zi which can annihilate efficiently 
into vector boson pairs, ameliorating the SUSY little hierarchy problem, while in accord 
with the measured abundance of cold dark matter in the Universe. While the analyses 
of 1^, ^ and take low M3 in an ad hoc fashion, the required gaugino mass pattern 
can also be obtained by allowing the SUSY breaking i^-term to transform as appropriate 
linear combinations of fields contained in the symmetric product of two adjoints of the 
unifying gauge group [p. We note here that a top-down scenario that naturally leads to 
low M3, low |^[ and light top squarks occurs in so-called mixed moduli-anomaly mediated 
SUSY breaking models, also referred to as mirage unification models, wherein moduli con- 
tributions give universal gaugino mass terms, but comparable gaugino mass splittings from 
anomaly-mediation reduce the value of M3, owing to the negative SU{3) beta function[^^. 

In this paper, we explore the phenomenological implications of compressed SUSY. We 
divide our discussion into two different model lines. In Case A (examined in Sec. ^), we 
adopt a model line from Ref. which is continuously connected to mSUGRA via variation 
of the gaugino mass M3, but with a non-zero Aq parameter. By dialing M3 to smaller 
values, the top squark mass is decreased, and the relic density is ultimately dominated 
by annihilation to tt via light ti exchange. The neutralino, however, remains essentially 
bino-like.^ The enhanced neutralino annihilation rate in turn implies an enhanced DM 
annihilation rate in the galactic halo|^, and we show that indirect DM search rates are 

'^For a further model with compressed spectra, see Bae et al, Ref. 

^If M3 is reduced farther, the neutralino develops a significant higgsino component and leads to mixed 
higgsino dark matter as already mentioned, unless of course, this range of A/3 is forbidden because ti 
becomes the LSP. 



1.5Mi = M2 = 3M3, 
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thus enhanced relative to mSUGRA. In addition, the low fi value and low rriq values typical 
of compressed SUSY result in enhanced rates for direct DM detection, and detection via 
muon telescopes. For this case, when the measured abundance of CDM is achieved, we 
arrive at a small mass gap solution where g — > tti dominantly, followed by ti — > cZi. 
In addition, the dominant decays Wi — > bii and Z2 Zih suggest that compressed 
SUSY LHC signatures are expected to be lepton poor, although robust rates for multi-jet 
+EJ^^^^ signals remain. We note, however, that Z2 — > ZZi has a branching fraction of a 
few percent. This, combined with the enormous rate for the production of sub-TeV scale 
gluinos (in the dark-matter-allowed regions) makes the multi-lepton signal observable in 
the cases we examined. 

In Case B (examined in Sec. ^), we consider a model line from MartinQ with l.bMi = 
M2 = 3M3. In this case as well, DM direct and indirect detection rates are larger than for 
the case of unified gaugino masses (with large |//|), and may possibly be detectable via ton 
size noble element detectors, or perhaps via anti-particle and gamma ray searches if the 
(currently undetermined) halo dark matter distribution turns out to be suitably clumpy, 
even though Zi remains dominantly bino-like. Since the mass gap m^^ ""^^i greater 
than rrih + Mw, we implement the 3-body decay ti bWZi into Isajet 7.76 (which we 
use for spectra and event generation). We find regions with a large branching fraction for 
tl bWZi decays, so that when this mode dominates, leptonic signals from gluino and 
squark cascade decays occur at observable levels. 



2. Case A: Low M3 scenario with continuous connection to mSUGRA 

In this section, we examine a model line based on mSUGRA, but with M3(GUT) as an 
independent parameter, with parameter space 

mo, mi/2, M3, Aq, tan/3, sign{ij,), (2.1) 

where we take the GUT scale values'^ Mi = M2 = mi/2 and adopt mj = 175 GeV to 
conform with Martin Q. The phenomenology of this scenario has been investigated in 
depth in Ref. for Aq = 0, where a low enough value of M3 <C mi/2 leads to a small 
fj, parameter, and hence the correct dark matter relic abundance via mixed higgsino DM. 
In the case studied here, we adopt a value of ^0 = —^■^^1/21 which helps reduce m^-^ 
compared with a choice of ^0 = 0, so that we can obtain dominant ZiZi annihilation 
into tt via a light ii exchange. Of course, if m^^ — m^^ becomes small enough, ii-Zi co- 
annihilation will also be important. Since for a bino-like LSP m^^ ~ 0.4mi/2, we will need 

mi/2 ~ 450 GeV so that m^^ > nit. Thus, we adopt mi/2 = 500 GeV, and take mo = 340 
GeV, tan/3 = 10 and ^ > in accord with Martin[0]. 



The mass spectrum - generated using Isajet 7.76|10] - is shown versus M3 in Fig. 
|l]a). In our illustration, M3 = 500 GeV corresponds to the mSUGRA model. Here, the 
spectrum shows the well-known feature that the colored sparticles (squarks and gluinos) 



^We will henceforth not explicitly specify the scale of the gaugino mass parameters, but this should be 
clear from the context whether we are referring to the parameters at the weak or at the GUT scale. 
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are split from, and much heavier than, the hghter uncolored sparticles. As M3 decreases 
from 500 GeV, the gluino, and via RGE effects also squark, masses drop giving rise to the 
"compressed SUSY" mass spectrum. The ti squark is the lightest of the squarks, owing 
to Yukawa coupling and intra-generational mixing effects, and its mass drops below 
and TTifj around M3 ~ 300 GeV. We note that the diminished squark masses feed into the 
Higgs soft masses via the RGEs, and give rise to a falling n parameter as M3 drops. The 
end of parameter space occurs at M3 ~ 238 GeV, where the ii becomes the LSP, and so is 
excluded by limits on stable charged or colored relics from the Big Bang. We see that not 
only m , but also , is significantly smaller than fi even at the lower end of M3 where 
the WMAP constraint is satisfied, so though Zi develops a significantly larger higgsino 
component compared to mSUGRA where it retains its bino-like character. 

In Fig. , we show the neutralino relic density h? versus M3 for the same param- 



eters as in frame a), using the IsaReD program 1 11]. For the mSUGRA case of M3 = 500 



GeV, ^ model would be cosmologically excluded, at least if we 

assume thermal relics and standard Big Bang cosmology. As M3 decreases from 500 GeV, 
il^^/i^ drops slowly until below M3 ~ 300 GeV a more rapid fall-off brings into 
accord with the WMAP measurement, which occurs for M3 ~ 255 GeV. At this point, the 
ti is rather light, with ~ 200 GeV, and m^^ ~ 230 GeV. 

In Fig. |2|, we show the integrated thermally weighted neutralino annihilation cross 
sections times relative velocity versus M3 as obtained using IsaReD, for various neutralino 
annihilation and co-annihilation processes. Here, x is the temperature in units of the 
LSP mass. The neutralino relic density is determined by the inverse of the sum shown 
by the solid red line, so that large annihilation cross sections yield low relic densities. In 
the mSUGRA case with M3 = 500 GeV, the neutralino annihilation rate is dominated 
by annihilation to leptons via t-channel slepton exchange. As M3 decreases, the squark 
masses, and especially the ti mass, decrease, so that ZiZi ti becomes increasingly 
important, and in fact dominates the annihilation rate for 240 GeV < M3 < 340 GeV. For 
lower M3 values, the ii — Zi mass gap is below 30 GeV, and top-squark co-annihilation then 
dominates, although in this narrow range ^^x^^ does not saturate the measured CDM relic 
density. We also see that as M3 decreases, annihilation to WW, ZZ and hh also increase 
in strength due to the lower /i value, and increasing higgsino component of the neutralino^. 
However, these channels never dominate in this case. 

In compressed SUSY, a light top squark is desirable in that it enhances the neutralino 
annihilation rate, and brings the relic density prediction into accord with observation, 
providing yet another mechanism for reconciliation of the predicted DM relic density with 
observation. However, generically a light top squark also enhances SUSY loop contributions 
to the decay b — > S7|T^. In Fig. |3|, we show the branching fraction BF{h — > 57) vs. M3 for 
the same parameters as in Fig. ||. In the mSUGRA case, the predicted branching fraction is 
in accord with the measured value: BF{b sj) = (3.55 ±0.26) x 10"^ from a combination 
of CLEO, Belle and BABAR data||l3|. However, the light ii in the low M3 region reduces 



*We have traced the turnover at low M3 in the various curves to a drop in the freeze out temperature 
that determines the range of integration. 
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the branching fraction well below the measured value. Of course, this branching fraction 
is also sensitive to other model parameters, e.g. tan/3. The point, however, is that for the 
light ii case, the SUSY contribution is generically comparable to the SM contribution, so 
that these must fortuitously combine to be consistent with the experimental value, which 
itself is in good agreement with the SM prediction. At the very least, in the absence of any 
real theory of flavor, (such fortuitous) agreement with the measured value, which agrees 
weh with the SM prediction [|14| BF{b sj) = (3.29±0.33) x 10"^, can always be obtained 
by allowing a small flavor violation in the soft parameter matrices at the GUT scale. 

2.1 Scenario A: dark matter searches 

Next, we investigate prospects for dark matter searches for the case A model line. We first 



calculate the spin- independent neutralino-proton scattering cross section using IsaReS|15|, 
and plot the results in Fig. ^a). In the case of mSUGRA at M3 = 500 GeV, the cross section 
(75/ (Zip) ~ 10"^'' pb, which is near the projected limit of future ton-scale noble liquid dark 
matter detectors. As M3 decreases, the squark masses also decrease, which increases the 
neutralino-proton scattering rate, which occurs primarily via squark exchange diagrams. 
Furthermore, a reduced value of |/u| is obtained for the low value of IM3I, resulting in an 
increased higgsino component of Zi (which still remains bino-like) so that the contribution 
to the direct detection cross section via the Higgs exchange diagram is correspondingly 
increased. By the time we reach ^^_^h? ~ 0.1 at M3 ~ 255 GeV, the direct detection cross 
section has grown by an order of magnitude, to just above 10~^ pb. This is a general 
feature of models with a low M3 value|p: for a given (bino-like) neutralino mass, direct 
detection rates are enhanced in the low M3 case. 

In Fig. ^6), we show the flux of muons expected to be measured at a neutrino telescope 
from neutralino annihilation into muon neutrinos in the core of the sun. In this and other 
indirect detection rates, we have implemented the Isajet/DarkSUSY interface []T6| . We re- 
quire muons to have energy > 50 GeV, the threshold for the IceCube detector[17|. 



In this case, the rate is again enhanced in going from mSUGRA to compressed SUSY, 
primarily because of the diminution of squark mass and the reduced value of as al- 
ready discussed above: these increase the spin-dependent neutralino-nucleon scattering 
cross section and enhance the IceCube rate because of the increased capture of neutralinos 
by the sun. In the WMAP-allowed region, the neutralinos mainly annihilate to tt pairs, 
so that the energy of the neutrino from top decays is shared with the accompanying b and 
the daughter muon. We see that although the flux of muon neutrinos corresponding to 
Efj_ > 50 GeV increases by a factor of ~ 500 in going from mSUGRA to the compressed 
SUSY case illustrated here, the flux of muon neutrinos is still below the reach of IceCube, 
primarily because the neutralino is still mostly bino-like. 

For positrons and anti-protons, we evaluate the averaged differential antiparticle flux 
in a projected energy bin centered at a kinetic energy of 20 GeV, where we expect optimal 



statistics and signal-to-background ratio at space-borne antiparticle detectors 1 18]. We take 



the experimental sensitivity to be that of the Pamela experiment after three years of data- 



taking as our benchmark] 19 1 . The expected fluxes depend on the (unknown) details of 



the neutralino distribution in our galactic halo. Here, we assume a scenario where baryon 
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infall causes progressive deepening of the gravitational potential well, and a clumpy halo 
distribution is obtained: the Adiabatically Contracted A^03 Halo Model [pO|] 

In Fig. m we show the expected positron flux in frame a) and the expected anti-proton 
flux in frame b) versus M3 for the same parameters as in Fig. ||. We see that in each 
case the antimatter flux jumps by a factor of 10^ in going from mSUGRA to compressed 
SUSY, largely due to the enhanced annihilation rate into tt pairs, and the concomitant hard 
spectrum of e"'"s and ps that ensue. In the case shown, the positron flux for compressed 
SUSY is somewhat below the Pamela reach, while the p flux is near the Pamela reach. We 
warn the reader that for the smooth Burkert halo profile [21| the signals are significantly 
smaller and beyond the sensitivity of Pamela. 

We have also evaluated the average differential anti-deuteron flux in the 0.1 < Tg < 
0.25 GeV range, where Tg stands for the antideuteron kinetic energy per nucleon, and 
compared it to the estimated sensitivity of GAPS for an ultra-long duration balloon-borne 
experiment 1 22]. We see in Fig. |6|a) that the antideuteron flux is again enhanced by a 
factor of ~ 10^ in going from mSUGRA to compressed SUSY, and in fact moves above the 
detectability limit of the GAPS experiment. For the Burkert halo profile, the estimated 
flux for the WMAP-allowed range of M3 is essentially at the edge of detectability. 

Indirect detection of neutralinos is also possible via the detection of high energy gamma 
raysp3[ produced by neutralino annihilation in the center of our Galaxyp^]. These will 
also be searched for by the GLAST collaboration ||2^ . We have evaluated expectations for 
the integrated continuum 7 ray flux above an = 1 GeV threshold versus M3 in Fig. |6|6). 
These projections are extremely sensitive to the assumed neutralino halo distribution, and 
drop by more than four orders of magnitude for the Burkert halo profile. This makes it 
difficult to make any definitive statement about the detectability of this signal (which could 
serve to map the halo profile rather than a diagnostic of the nature of the DM particle). 
However, once again we see a factor of ~ 100 enhancement in detection rate in moving from 
the mSUGRA case where M3 = 500 GeV to the compressed SUSY case with M3 ~ 255 
GeV. 



2.2 Scenario A: LHC searches 

As Martin notes the compressed SUSY mass spectra are generally too heavy for suc- 
cessful sparticle searches at the Fermilab Tevatron. However, ( |1.1D implies an upper bound 
on the bino mass, and since we reduce M3 from its unified value, implies that gluinos must 
be relatively light so that multi-jet -|- multilepton +E^^^^ events from SUSY should be 
produced in abundance at the CERN LHC, due to turn on in 2008. In this section, we 
investigate the collider signals expected after cuts for various signal topologies at the LHC. 

At the CERN LHC, gluino and squark pair production will be the dominant SUSY 
production reactions. Gluino and squark production will be followed by their cascade 
decays |p6[] , resulting in a variety of events with jets, isolated leptons and missing energy. 
A large number of signals emerge, and can be classified by the number of isolated leptons 
present. The signal channels we examine include i.) no isolated leptons plus jets plus E^^^^ 
{0£), a.) single isolated lepton plus jets plus E'™'^'^ (1£), in.) two opposite sign isolated 
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leptons plus jets plus E^^^^ (OS), iv.) two same sign isolated leptons plus jets plus E^^^^ 
(SS) and v.) three isolated leptons plus jets plus E^^^^ {3i). 

The reach of the CERN LHC for SUSY has been estimated for the mSUGRA model 
in Ref. p^, 28| for low values of tan/3 and in Ref. |29| for large tan/3 values. We adopt the 



cuts and background levels presented in Ref. |27| for our analysis of the signal channels 
listed above. Hadronic clusters with Et > 100 GeV and |?7(jet)| < 3 within a cone of 
size R = + A02 = 0.7 are classified as jets. Muons and electrons are classified as 

isolated if they have pT > 10 GeV, \r]{i)\ < 2.5, and the visible activity within a cone of 
R = 0.3 about the lepton direction is less than £'r(cone) = 5 GeV. 

Following Ref. [^], we required that the jet multiplicity, njet > 2, transverse sphericity 
St > 0.2, Eriji), and further, that £^t(J2) > and EJ^''''' > Elf,, where the cut 

parameter Elp is chosen to roughly optimize the signal from gluino and squark production. 
For the leptons we require pt{£) > 20 GeV {£ = e ot h) and Mt{£, EJ^''''') > 100 GeV for 
the li signal. For the OS, SS and 3i channels, we require that the two hardest leptons 
have pt > 20 GeV. We have also applied a cut on the transverse plane angle A(j){E^^^^ , jc) 
between and closest jet: 30° < Acf) < 90°, in the case of the E'™'^'^ channel, i). 

Our results are shown in Fig. for a rather loose choice of the cut parameter E^ = 
100 GeV. We see that as M3 decreases from the mSUGRA value of 500 GeV, the signal 
cross sections increase. The increase is mainly due to increased total gluino and squark 
production cross sections, due to their decreasing masses. When we reach the DM -allowed 
compressed SUSY spectrum at M3 ~ 250 GeV, however, the leptonic signals suffer a steep 
drop-off, while the E^^^^+ jets signal increases somewhat. This is due to the fact that in 
this case, Wi — > bti turns on and dominates the Wi branching fraction, while ti cZ\ 
at essentially 100%. Thus, no isolated leptons come from chargino decay. Likewise, Z2 — > 
Z\h at around 90% branching fraction, so isolated leptons from Z2 decays come from the 
subdominant decay chain Z^ — > Z\Z which has a branching fraction of ^ 8%. Isolated 
leptons still arise from g — > ti\ decay, followed by semi-leptonic top decay, but in general, 
we expect in compressed SUSY models with a small t\ — Z\ mass gap and fn^^ > m^^ +mb 
that the fraction of signal events containing isolated leptons will be much lower than the 
usual prediction from models like mSUGRA with gaugino mass unification. We regard a 
signal to be observable if for the given integrated luminosity, i) the statistical significance 
of the signal exceeds 5c7, ii) S/B > 0.25, and iii) S > 10 events. The minimum observable 
cross sections for each topology are shown by the dashed horizontal bars in the figure. We 
see that even for the low value of E^ = 100 GeV, all but the opposite sign dilepton signal 
should be observable with an integrated luminosity of 100 fb^^, and frequently even with 
a much lower integrated luminosity, at least for parameters in the WMAP-allowed region. 
Although we do not show this, we have checked that with E^ = 200 GeV, the OS signal is 
easily observable,^ and furthermore, the Qi signal is not as close to the observability limit. 



^Since the OS dileptons come primarily from the decay of an on-sheU Z boson, it is possible that this 
signal would actually be observable even for — 100 GeV. 
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3. Case B: Non-universal gaugino masses and a large mass gap 



In this section, we explore Case B, the compressed SUSY model line originally suggested 
by Martin where at Q = Mgut, 1-5Mi = M2 = 3M3, with mo = 340 GeV, ^0 = -0.75Mi, 
tan /3 = 10 and /i > 0. We first display the variation of the sparticle mass spectrum with 
Ml in Fig. ^a). The upper end of parameter space is limited by Mi ~ 1000 GeV, where for 
higher Mi values the ti becomes the LSP. This implies an upper bound of 1200 GeV (1100- 
1400 GeV) on gluino (squark) masses, ensuring their copious production at the LHC. The 
lower range of Ml is bounded by Mi ~ 160 GeV, since for lower Mi values, the value of ra~ 
drops below limits from LEP2 [^]. In the intermediate region with 440 GeV < Mi < 1000 
GeV, the ti is relatively light, and is the next-to-lightest SUSY particle (NLSP). More 
importantly from our perspective, in this mass regime > nit, ensuring that ZiZi ti 
was accessible in the early Universe. 

In Fig. ^6), we show the neutralino relic density as a function of Mi for the same 
parameters as in frame a). There is a wide range of Mi : 400 — 800 GeV where the relic 
abundance is in close accord with the WMAP measured value. It should be possible to 
bring this in accord with the WMAP value by slightly tweaking ^o- For yet larger values 
of Ml, tiZi and titi annihilation rates become large, and the Zi relic density no longer 
saturates the observed density of CDM; i.e. the DM would be multi-component in this 
case. In contrast, when Mi drops below ~ 400 GeV, corresponding to m^^ < m^, the 
prediction for ^^Ji^ rises above the WMAP measurement, excluding Zi as a thermal relic. 
For Ml ~ 150 GeV ~ a range excluded by the LEP2 chargino mass limit - there is a 
double dip structure where 2m^^ ~ rrih or M^, and so neutralinos can efficiently annihilate 
through these s— channel poles. 

In Fig. ^, we show the integrated thermally weighted neutralino annihilation (and 
co-annihilation) cross section times relative velocity versus Mi for the same parameters as 
in Fig. 1^. For Mi ~ 750 GeV, the ti — Zi mass gap is so low that Ziti co-annihilation, 
and eventually titi annihilation (not shown in the figure), dominates and we get too small 
a relic abundance. In the range 400 GeV ~ Ml ~ 750 GeV, Zi^i ^ tt dominates, 
so agreement with the relic density is obtained as envisioned by Martin ||2|. For Mi ~ 
400 GeV, annihilation into ti is not allowed (except for Zis in the high energy tail of 
the thermal distribution), and so annihilation takes place dominantly into WW (via the 
higgsino component) and into quarks and leptons. At the h and Z poles (inside the LEP- 
forbidden region), annihilation into down- type fermions dominates. 



The branching fraction BF{h — > 57) is shown in Fig. 10 versus Mi for the same 
parameters as in Fig. ^. Here we note that for large Mi, the branching fraction is close to- 
albeit somewhat below- its measured value. However, as Mi decreases, the ti and Wi both 
become lighter, and SUSY loop contributions to the branching fraction move the predicted 
branching fraction away from its observed value. In this case, as in Sec. we would expect 
a somewhat suppressed value of BF{h 57) compared to its SM predicted rate. We recall 
as before that it should be possible to bring this result into accord with experiment by 
allowing, for instance, some flavor structure in the soft SUSY breaking sector. 
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3.1 Case B: dark matter searches 

Here, we examine direct and indirect dark matter detection rates for the compressed SUSY 
model hne Case B. We begin by considering the prospects for direct detection in Fig. |ll|a) 
where we show the spin-independent Z\-p cross section as a function of the GUT scale 
parameter M\. The cross section increases as M\ decreases due to decreasing squark 
masses, and a decreasing value of the /U parameter. The range relevant for compressed 
SUSY, Afi : 400 - 750 GeV, has asi[ZxV) ~ 5 - 15 x 10"^ pb, which at its high end 
is within an order of magnitude of the current limit from XENON- 10 [^], and should be 
detectable by SuperCDMS or 100-1000 kg noble liquid DM detectors. Projections for direct 
detection are somewhat more optimistic than in Case A, mostly because the value of /u is 
relatively smaller in Case B. 

In frame 6), we show the flux of muons with > 50 GeV expected at neutrino 
telescopes due to neutralino annihilation in the solar core. As Mi decreases from 1000 
GeV, the rate slightly increases, due to an increasing spin-dependent Zi-nucleon scattering 
rate, but for the Afi range of interest, the flux remains somewhat below the IceCube 
detectable level. For Mi < 400 GeV, the rate jumps to higher levels. This jump can 
be understood from Fig. from which we infer that since the neutralino capture and 
annihilation processes are in equilibrium, the fraction of captured neutralinos that directly 
annihilate into z^P jumps once annihilation to tt turns off, and it is these very high energy 
neutrinos which have the greatest chance of being detected by IceCube. For Mi > 400 
GeV, ZiZi annihilates mainly into tt, and the fraction of direct neutralino annihilation 
into neutrinos is lower. 

In Fig. H we show the flux of a) positrons and b) anti-protons from neutralino an- 
nihilations in the galactic halo expected in Case B versus Mi, assuming the clumpy halo 
as given by the adiabatically contracted N03 halo model. We evaluate the signal in the 
same energy bins and apply the same sensitivity criterion as in Fig. ^. The flux of e^s is 
everywhere below the Pamela sensitivity even for our favorable choice of halo distribution. 
However, the results do show some structure and enhancement in the compressed SUSY 
range of Mi : 440 — 800 GeV. In this regime, > mt so that ZiZi — > tt can occur in 
the present galactic halo as well as in the early Universe. The turn-on of the tt annihila- 
tion mode is clearly seen at Mi ~ 440 GeV. In the case of the p flux, the signal actually 
increases enough to suggest some range of observability at Pamela. 

In Fig. H we show a) the flux of anti-deuterons along with the reach of the GAPS 
experiment, and b) the flux of gamma rays from the galactic center with E^^ > 1 GeV. In the 
case of Ds, the entire compressed SUSY range is above the GAPS sensitivity. We caution, 
however, that for the smooth Burkert halo profile, projections could be lower by a factor 
10-15. For 7s, the entire range should be probed by GLAST, although these projections 
are extremely sensitive to our assumed halo distribution; for other halo choices - such as 
the Burkert profile, the rates scale downwards by over four orders of magnitude, and could 
fall below the projected sensitivity of GLAST. However, in both the case of Ds and 7s, 
a sharp increase in indirect detection rate occurs when ZiZi tt turns on at Mi ~ 440 
GeV. 
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3.2 Scenario B: LHC searches 

An important issue for evaluating collider signals in compressed SUSY with a large ti — Z\ 
mass gap is to evaluate the t\ decay branching fractions correctly when the t\ is the NLSP. 
In this case, the flavor changing decay t\ — > cZ\ may compete with the three-body decay 
t\ — > bWZi if the latter decay mode is kinematically allowed. We implement the three 
body decay into Isajet 7.76 using the squared matrix element calculated by Porod and 
Wohrmann|32|.^ We also update the Isajet formulae for the flavor changing two-body 



decay using the one-step integration approximation of Hikasa and Kobayaship^], but with 
the correct neutralino eigenvectors p5[ |. We have checked that this single step integration 
systematically over-estimates the width of the loop decay ti cZi, so that if we find the 
three body decay to be dominant within our approximation, this will be the case also with 
the correct calculation. In Fig. 0, we show the branching fraction of ii versus Mi for the 
same parameters as in Fig. ^ We see that at large Mi where m^^ < rrih + Mw + ni^^ , the 
ii decays entirely into cZi^. For lower Mi values, the ii — > bWZi decay mode opens up 
and in fact dominates the two-body mode for Mi : 400 - 460 GeV. For Mi ~ 400 GeV, 
then m^^ > rrif, + tn^^ , so that ti bWi turns on and dominates the branching fraction. 
In this regime, for this case at least, < mt, so this range is not as interesting from 
the perspective of obtaining agreement with the relic density measurement via neutralino 
annihilation to top quarks. 

Once the correct decay patterns are implemented, we can generate collider events and 
examine signal rates after cuts. We present multi-lepton plus multi-jet cross sections 



at the LHC for Case B, using the same analysis as in Sec. 2^, except with = 200 GeV. 



The results are shown in Fig. 15. For low Mi values, the squark and gluino masses are 
relatively light, and SUSY particle production cross sections are large at the CERN LHC. 
Nevertheless, signals in the Oi, li and OS channels fall below the S/B > 0.25 level for 
Ml in the interesting range of 400-800 GeV if we choose = 100 GeV. In contrast, with 
E"^ = 200 GeV shown in the figure, signals in all channels are observable for the entire 
range of Mi.*^ As Mi increases, the g — Wi and q — Wi mass gaps actually increase, and 
we get an increase in the multi-lepton signal rates. These show a rapid drop off beyond 
Ml = 450-500 GeV, where ti produced in gluino cascades decay via ti — > cZi. There is 
no analogous drop-off in the Oi or even in the li channels (since it is not difficult to get 
a single lepton somewhere in the cascade, e.g via the decay of t). Thus, once the stop is 
light enough so it can only decay via ti — > cZi (which is the case over most of the Mi 
range of interest), we see a relative reduction of multi-leptonic signals compared with those 
containing just -E'^^'^'^-l-jets. Distinguishing Case B (with ti cZi) from Case A will be 
challenging at the LHC, but should be straightforward at a TeV linear collider. 



®We have, however, made one correction from their erratum. The term 2m~ \2pi, ■ pw + "i- j in 

Eq. (A. 2) should be replaced by 2m|^pi, • pw and not by 4m^^pi, ■ pw as stated in their erratum. 
M. MiihUeitner (private communication) has independently confirmed this factor, which also appears cor- 
rectly in the program SDECAY[|3|. 

^There are four body decay modes such as ii bff'Zi (where / is a SM fermion) which we have not 
evaluated, but which we expect to be smaller than the two-body decay 

*The — 100 GeV is better optimized for the signals in the SS and 3£ channels. 
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4. Summary and conclusions 



The generic prediction of the neutrahno rehc density from SUSY models falls somewhat 
above the measured value if sparticles are significantly heavier than ~ 100 GeV, as is 
likely to be the case given the direct constraints from LEP2 and the Tevatron, and in- 
direct constraints from low energy measurements. Within a particular framework such 
as mSUGRA, this means special regions of parameter space where at least one of: neu- 
tralino co-annihilation with staus/stops, neutralino resonance annihilation via A/H (this 
requires large tan /? values) or even h, or mixed higgsino DM in the hyperbolic branch/focus 
point region at large mo, obtains. Each of these alternatives would have implications for 
SUSY signals both at colliders, as well as for direct and indirect searches for DM. Unfor- 
tunately, these implications are not robust to small changes in the model. Allowing for 
non-universality of gaugino or Higgs scalar mass parameters leads to one-parameter exten- 
sions of the mSUGRA model where the implications of the WMAP measurement can be 
strikingly different. For instance, non- universal Higgs mass models allow mixed Higgsino 
DM for low values of rriQ, and Higgs resonance annihilation for all values of tan /? [37|. Non- 
universal gaugino masses allow new possibilities, such as mixed wino DM or bino-wino 
co-annihilation ||3^ that are precluded in models with unified gaugino masses (but realized 
in other frameworks). These studies suggest that it would be premature to blindly use the 
measured relic density to make definitive projections for what should/should not be seen 
at the LHC or in DM searches. Already there exist numerous alternatives (with different 
phenomenological outcomes) to choose from, and only experiment can zero in on nature's 
choice. 

In this vein, Martin Q recently pointed out yet another possibility to obtain agreement 
with the observed CDM relic density. He noted that if m^^ > rrit and ii is not much heavier 
than Zi, then ZiZi — > tt mediated by ti exchange (this process does not suffer the large 
p-wave suppression on account of the large top mass) may dominate in the early universe. 
This scenario can be realized for a bino-like Zi only if gluinos (and through the RGEs, also 
the squarks) are not very heavy, leading to a "compressed SUSY" spectrum. In this paper, 
we have examined two different model lines that realize Martin's idea, and quantified the 
implications for SUSY searches at the LHC as well as via direct and indirect searches for 
DM. 

The first model line that we refer to as Case A is continuously connected to mSUGRA, 
and is in a sense an extension of our earlier work that we referred to as low IM3I dark 
matter, where relaxing the gaugino mass unification condition and allowing [Ms (GUT) | to 
be smaller than Mi ~ M2 led to viable solutions with mixed higgsino DM [||. In these 
studies, we used ^0 = for simplicity. Here, we choose instead ^0 = —1 -5^1/2, and 
lower M3 as before. This choice of Aq leads to a reduction in m^^, and remarkably, as 
M3 is reduced, m^^ becomes close to m^^, so that ZiZi — > tt can indeed be the dominant 
mechanism in the early Universe, with Zi retaining its bino-like character. While the 
reduced gluino, and concomitantly squark, masses and the /i parameter, imply larger direct 
and indirect detection rates vis a vis models with gaugino mass unification, these rates are 
not large, primarily because the neutralino remains bino-like. Nevertheless ton size noble 
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liquid detectors should be able to directly see a WIMP signal (at least for parameters 
that give the observed relic density), while indirect searches for anti-deuteron at GAPS 
or gamma rays from our galactic center by GLAST may yield an observable signal, but 
only if the DM is favorably clumped. We project that there will be no detectable signal in 
Pamela or in IceCube. The scenario implies that gluinos and squarks cannot be too heavy 
so that the LHC should be awash in SUSY events, and the signal should be extricable from 
SM backgrounds with simple cuts. The characteristic feature of the scenario is the relative 
reduction of the signal in multi-lepton channels relative to that in Oi or li channels. The 
large production rate nevertheless implies there should be an observable signal in all the 
channels shown in Fig. 0. A significant fraction of OS dilepton and trilepton events may 
contain a real Z boson. 

The second model line that we examine (and refer to as Case B) is the one suggested 
by Martin in his original proposal. Here, we adopt non-universal boundary conditions 
l.bMi = M2 = 3M3 for the GUT scale gaugino mass parameters. Prospects for direct 
detection may be somewhat better in this scenario: in favorable cases, the signal cross 
section may be just an order of magnitude away from the current upper bound. Indirect 
detection prospects are similar to those in Case A. There is no detectable signal at IceCube, 
potentially observable signals in GLAST or GAPS for favorable halo distributions, and 
possibly a marginal signal from p in Pamela. Experiments at the LHC should be able 
to detect a signal in all channels, albeit with somewhat harder cuts than in Case A, as 



illustrated in Fig. |15|. As in Case A, over most of the parameter range compatible with the 
relic density measurement, multi-lepton signals will occur at smaller rates. 

A light ti is the hallmark of the scenario. While its direct detection is not easy at 
the LHC,^ its presence along with that of a not-too-heavy chargino leads to a significant 
SUSY contribution to the 6 — > 57 branching ratio, and likely also to the branching ratio 
and distributions for b — > sii decays (that we have not examined). Indeed, for both cases 
that we examined, the former turns out to be smaller than its measured value. While it is 
certainly true that we can always reproduce the observed branching fraction by tweaking 
the flavour structure of soft-SUSY-breaking parameters, it would seem unlikely this would 
be "just right" to yield the SM prediction. It, therefore, seems that a deviation of the 
patterns of rare flavor- violating decays of 6-quarks from SM expectations should generically 
be expected in these scenarios. 
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mg=340 GeV, mj^2=500 GeV, -l.5m^^^, tan(3=10, \i>0, m=175 GeV 

10" I IM I I I I 11 




M3 (GeV) 

Figure 1: a) Sparticle mass spectrum for the case with toq = 340 GeV, Mi = M2 = 500 GeV, 
Aq = — 1.5toi/2; tan/3 = 10, /i > and nit — 175 GeV, versus GUT scale SU{3) gaugino mass 
parameter M3, and b) neutraUno rehc density versus for same parameters as in frame a). 
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Figure 2: Integrated thermally weighted cross sections times relative velocity for processes that 
may be relevant for the calculation of the Z\ relic density in the Big Bang versus M3. We illustrate 
these for the same parameters as in Fig. 
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mj=340 GeV, m|,,=500 GeV, A|= -LSm,,,, tanP=10, n>0, m=175 GeV 



m,=340 GeV, m„, =500 GeV, A, = -1.5ra„,, tanP=10, n>0, m =175 GeV 
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Figure 4: a) Spin-independent neutralino-proton scattering cross section and h) flux of muons 
with i?^ > 50 GeV at IceCube versus Afs for same parameters as in Fig. |l|. 
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m,=340 GeV, m^,, =500 GeV, Aj= -1.5m„„ lanP=10, n>0, m=175 GeV 
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Figure 5: a) Expected positron flux and h) antiproton flux versus M3 for same parameters as in 
Fig. 0. The dashed line shows the expected three year sensitivity of Pamela. 



mj=340 GeV, m|,,=500GeV, A^= -1.5m|,„ tanP=10, n>0, m=175 GeV 



mj=340 GeV, m,,, =500 GeV, A„= -1.5m|,„ lanP=10, n>0, m =175 GeV 
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Figure 6: a) Expected anti-deuteron flux and h) gamma ray flux versus A/3 for same parameters 
as in Fig. |^. The horizontal lines show the projected sensitivities of the GAPS and GLAST 
experiments. 
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mQ=340 GeV, m^/2=^00 GeV, A =-150 GeV, tanp=10, ^i>0, m = 175 GeV 
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Figure 7: Signal rates at the CERN LHC for various multi-jet plus multi-lepton -\-E^^^^ event 
topologies after cuts listed in the text with the cut parameter = 100 GeV versus M3 for same 
parameters as in Fig. 0. The horizontal dotted lines show the minimum observable cross section 
for E"^ ~ 100 GeV, assuming an integrated luminosity of 100 fb^^. 
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mg=340 GeV, A„= -0.75Mj, tan(5=10, n>0, m=175 GeV 
1.5Mj = M, = 3M, 
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mQ=340 GeV, Aq=-0.75Mj, tan|3=10, }i>0, m =175 GeV 
1.5M| =M2 = 3M3 




Mj (GeV) 

Figure 8: a): Sparticle mass spectrum as a function of the GUT scale gaugino mass Mi for 
Case B, where mo = 340, 1.5Mi = Ma ^ SAh GeV, Aq = -0.75Mi, tan/3 = 10, ^ > and 
mt = 175 GeV. b): Neutralino relic density versus AIi for same parameters as in frame a). 
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mo=340 GeV, Ao=-0.75M^, tanp=10, ^i>0, m,=175 GeV:1 .5M^=M2=3M3 
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Figure 9: Integrated thermally weighted neutralino annihilation (or co-annihilation) cross sections 
times relative velocity, for same parameters as in Fig. ^ versus Mi. The processes shown do not 
saturate the total at very large values of Mi because we have not plotted titi annihilation which 
becomes very important there because ti becomes very close to m^^ at the upper end of the Mi 
range. 



m^=3W GeV, A^= -0.75Mj, tanP=10, H>0, m =175 GeV 
1.5Mj =M2 = 3M3 




Figure 10: Branching fraction for b ^ decay versus Mi for Case B, for the same model 
parameters as in Fig. |^. 
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Ti„=340 GeV, A|j=-0.75M|, tanp=10, Jl>0, m =175 GeV 
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m|,=340 GcV, A|j=-0.75M|, timp=10, n>0. ra =175 GcV 
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Figure 11: a) Spin- independent neutralino-proton scattering cross section and b) flux of muons 
with Efj^ > 50 GeV at IceCube versus Mi for Case B with the same parameters as in Fig. 
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Figure 12: a) Expected positron flux and b) antiproton flux for Case B versus Mi for same 
parameters as in Fig. R 
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Figure 13: a) Expected anti-deuteron flux and 6) gamma ray flux for the compressed SUSY 
Case B versus Mi for same parameters as in Fig. @. 
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mg=340 GeV, Aq=-0.75Mj, tanP=10, \l>0, m =175 GeV 
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Figure 14: Branching fraction of ii vs. GUT scale gaugino mass Mi for same parameters as in 
Fig. |. 



mQ=340 GeV, Aq=-0.75Mj, tanp=10, ^>0, m = 175 GeV 



10 



10 



10 



10^ 



^ 10' 



10 



10 



10" 



Observability 
wvith 100 fb' 

,..0l 

'"11 



• OS 



.ss 

■31 



l.SMj =M2 = 3M3 



-p-i— 1— 1— 1— p-i— 1— 1— 1— p-i— 1— 1— 1— p 



-p-i— 1— 1— 1— i- 



E\= 200 GeV 




100 200 300 400 500 600 

Mj (GeV) 



700 



800 



900 



Figure 15: Signal rates for Case B for various multi-jet plus multi-lepton +E^^^^ events at the 
CERN LHC, after cuts detailed in the text and = 200 GeV, versus Mi for same parameters as in 
Fig. ||. The horizontal dotted lines show the minimum observable cross section for i?^ = 200 GeV, 
assuming an integrated luminosity of 100 fb^^. 
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